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Abstract 

It has been suggested that the observed pulsar velocities are caused by an asym- 
metric neutrino emission from a hot neutron star during the first seconds after the su- 
pernova collapse. We calculate the magnitude of gravitational waves produced by the 
asymmetries in the emission of neutrinos. The resulting periodic gravitational waves 
may be detectable by LIGO and LISA in the event of a nearby supernova explosion. 

1 Introduction 

It has been suggested that the proper motions of pulsars |lj may be caused by a slight 
assymetry in the ejection of neutrinos |21 El E], which carry away most of the energy of 
a supernova. In fact, it would take only a 1 % assymetry in the neutrino emmission to 
achieve the observed neutron star velocities. Such an assymetry could be caused by neutrino 
oscillations in matter, in a strong magnetic field. Although neutrino magnetic moments 
are negligible, the neutrino interactions with matter depend on the polarization of matter 
fermions The viability of the neutrino kick mechanism has been discussed at length in 
several papers |21 El El- This mechanism can work if there is a sterile neutrino with mass 
1-20 keV and a small, sin 6 ~ 10~^, mixing with the electron neutrino. It is intriguing, that 
the same neutrino could make up the cosmological dark matter jl] . Unfortunately, detection 
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of such a neutrino in a terrestrial experiment is probably not feasible at present because the 
mixing angle sin 9 is very small. However, in the event of a nearby supernova, the asymmetric 
emission of neutrinos from a rotating nascent neutron star can produce gravity waves that 
may be detectable. 

In this paper we will describe the expected gravity waves (GW) from a pulsar kick due to 
neutrino conversions taking into account the rotation of the neutron star and the neutrino 
asymmetry associated with this mechanism. Detection of such a signal can help understand 
the origin of the pulsar kicks. 

Gravitational waves from the supernova neutrinos have been considered in the literature jHl 
13 IH]. However, the source geometry was assumed somewhat different from that specific to 
the neutrino kick mechanism. Epstein jB] has calculated the burst radiation assuming the 
distribution of neutrinos is roughly elhptical. The assymetry was assumed to arise due to 
flattening by centripetal motion. Burrows and Hayes [7] performed a hydrodynamic simu- 
lation of a supernova with some imposed initial assymetry. Cuesta |8j described the burst 
radiation from the sudden neutrino conversions. He has also proposed that the GWs them- 
selves may be the source of the neutron star momentum kick, through their radiation back 
reaction 0. 

However, there is an additional, and possibly stronger, source of gravitational radiation 
due to the asymmetric emission of neutrinos from a rapidly rotating neutron star. One can 
think of the escaping neutrinos as forming an off-centered rotating beam, at some angle to 
the rotation axis, superimposed on a spherically symmetric distribution. The symmetric 
component does not produce gravity waves, but the 1% of all neutrinos that make up the 
asymmetric "beam" are a source of a periodic GW signal. 

2 Gravity waves from a spinning neutrino beam 

We will adopt Epstein's notation |6j, with a few minor modifications. We define the nontensor 
potentials 



where g^^, is the metric. Indeces of all non-tensors are raised and lowered with the Minkowski 
metric 77°"^ = rjap = diag(— 1, -t-1, Also, we define ,y = d/dx" , and take G = c = 1. 
This allows us to write the Einstein field equations as 






(2) 
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with the gauge condition 

e^\, = Q. (4) 

Here we followed Epstein in defining 

r^^^ = r'^f^ = T'^^ + 1^" (5) 

T^^^ is the matter stress-energy tensor and t^^ is the suitably defined gravitational stress- 
energy pseudo-tensor. 

Equation Q is in the form of a wave equation in fiat space. Its solution for outgoing 
wave boundary conditions is 

^^'^(t,x)=4 y"r'^^(t- |x-x'|,x')|x-x'|-MV (g) 

at any field point = (t,x) = (x", x*), i, j, ... = 1,2,3. We will assume all gravitational 
fields are weak enough that we may evaluate equations © and (jH} to leading order in 9^'^. 
Following Epstein [H], let us consider a source of the form 

r^^(t,x) =n'n^r-^cx{t-r)f{Q,t-r), (7) 

where n = x/r, r = |x|, f{Q,t) > 0, / f{Q,t)dQ = 1. This form describes neutrinos 
continuously released from x = at the speed of light. a(t) and f{fl,t) are the rate of 
energy loss (luminosity) and angular distribution, respectively. In the weak field limit, it 
is safe to treat all neutrinos as orriginating from x = and to use t — r =const as an 
approximation to a null geodesic. 

To obtain the detectable radiation we must find the transverse traceless (TT) part of the 
source. This can be done either with a projector or with a gauge transformation. Let us use 
the former. The detectable radiation in then 

^^^(t,x)=4 yr^^(t- |x-x'|,x')|x-x'rMV, (8) 

where 

r^T = {PlPi - \P''PuiV' (9) 



with 

= 5] - fi'fij, fi = . (lo) 



X — x' 



Here n is a unit vector in the x direction, n' is in the x' direction and n is in the x — x' 
direction. Let 9 be the angle between n and n', is the angle between n and n and if) + 9 
is the angle between n' and n. 
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Then assuming x is along the z-direction, and n' has polar coordinates {9, 0), we find 

1 



{nV)TT 



{n'yn'y)TT 



sin^(?/' + 6) cos^ 0(1 + cos^ip) - 1 
sin^{%l) + 9) [sin^ 0(1 + cos'^%1)) - 1 



(n"'n'y)TT = - sin^(?/' + 61) sin0cos0(l + cos^ V), 
{n'^n"')TT = ^sm^i/jsm^{ip + 9), 



{n"'n'^)TT 



[n''n'y)TT 



sin ip cos cos sin^ {ip + 9), 
sin -0 cos -0 sin sin^ (-0 + 6') . 



Using the identities 



sin0^ sin 6* sin(-0 + 6*) 



X — X 



X — X' 



+ r'^ — 2rr' cos 6', 



we write equation (0) as 

/oo 
f{n\t')a{t')6{t-t' -r)dt'. 
-oo 

Inserting this into equation (jHl) we obtain 

{n'n^)TTf{Sl'.t')a{t') 



9'Mt,^) =4 



— oo Jin Jo 



X — X' 



We can now perform the r' integration using equations (fTHj) and (fT^ to yield 
^^r(t,x) =4 

Let us define 



oo Jin t — t' — r COS 6*' 



a 



1. 



Then we can write all the angles in terms of 9' and a 



sin -0 



X — X 



■ sin 9' 



+ 2a 



a2 + 2a(l - cos 9') + 2(1 - cos 9'\ 
4 



sin 6*', 



(11) 
(12) 

(13) 

(14) 

(15) 

(16) 
(17) 

(18) 
(19) 

(20) 



^{t-\-^-^\-t' -r')dr'd£)!di! . (21) 



(22) 
(23) 



■ f I I a'\ ^ -a' 2{a + l - cos9') . 

sm(w + ) = r sm = — -, — -, — sm . 

' |x-x'| a2^2a(l-cos^') + 2(l-cos0') 

We can then write 

Ji-K a + 1 — cos 9' 

Next we write equation as 

^^^(t,x) = 4r-i Q'^{a,t')a{t')dt'. (25) 

Let us note that all of the dependence on t (not t') is either in the integration limit, or in a. 
Thus, we can write the time rate of change of 6^^ as 

e'it, x) = 4r-^Q'^{0, t - r)a{t - r) + Ar'^ / ^ ^ ' ^ a{t')dt' . (26) 

J-oo da 

We see that there are two contributions. The first is due to the new neutrinos being 
ejected at t — r, the ones whose effect is just becoming visible. The other, is due to the 
change in the effect of previously ejected neutrinos, as their distance and angle change. 
However, this second term is suppressed by an extra factor of r, so it is negligible as long as 
the radiation lasts for a time that is short compared to the distance between the source and 
the observer. (In a supernova, the event is about 10 s and the distance is thousands of light 
years so we can safely ignore this term.) 

For the first term, since a = and ip = 0, the angular terms simplify greatly. All of the 
6^rp terms drop out and we find 

(^^^0)TT(x,t) = 

-(^^^o)TT(x,t) = 2r-V(t-r) / (1 + cos^') cos207(fi',t - r)rff]', (27) 
(r^o)rr = 2r-^a{t-r) f {1 + cose') sin 2(1)' f{Q' ,t - r)dQ' . (28) 



3 (5-function Distribution 

At this point we simply need to choose a distribution f{Q', t) to find the form of the radiation. 
The simplest is just a S function in a particular direction. In essence, we are modeling the 
slight neutrino excess to one side as a 5-function jet. In reality the 6 function introduces 
some spurious results, but it also provides a basis for more accurate results. Let 

f{n',t) = 6{cose{t) ~ cose'{t))6{(i){t) - (f)'). (29) 
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Then we find 



e-,OTT(x,t) = 

-^^^OTT(x,^) = 2r-^a{t-r){l + cose{t-r))cos2(j){t-r), (30) 
e^^y,^ = 2r-^a(t-r){l + cose(t-r))sm2(l)(t-r). (31) 

The most striking feature is that while the radiation is suppressed for ^ = tt when the 
neutrino jet is away from the the observer as expected, there is no suppression near ^ = 
when the jet is toward the observer. We will see later that this is one of the spurious effects 
of the 5-function. 

The neutrino jet's motion could be quite complex. We will focus on one simplified 
possibility where the jet precesscs at an angular frequency a; around an axis that is inclined 
at an angle a relative to the z-axis, and always makes an angle of x relative to that axis. For 
simplicity we will always assume that the axis of precession is in the x — ^-plane. Rotating 
it away from this plane can be handled most easily by rotating the observer's axes. 

If q; = then the radiation is simply 

r^,0Tr(x,t) = 2r-V(t-r)(l + cosx)cos2o;(t-r), (32) 
^^^,0Tr(x,i) = 2r-^a{t-r){l + cosx)sm2u{t-r). (33) 

Which is just circularly polarized gravitational radiation at an angular frequency of 2a; which 
is exactly what we would expect. 

If we move a away from the 2;-axis. Things get more complicated. By simply rotating 
the axes we can find 

cos^ = cosacosx — sinasinxcos^o, (34) 

tan0 = sn^x^ 

sin a cos x + cos a sin x cos i^o 

o , 1 - tan^ (f) 

cos 20 = ■ — 5—, (36) 

1 + tanV 
. 2tan0 

1 + tan"^ 

where we have let 0o = (^{t—r). Using these gives us gravitational radiation that is elliptically 
polarized and and contains components at all angular frequencies nu where n is an integer. 
The radiation is generally dominated by lu and 2a; components, and in fact we will find that 
the higher frequency terms are spurious effects of the 5-function. 
The magnitude of any particular Fourier component is then 
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L is the luminosity of neutrinos, e is the excess in one direction (about 1%), oj is the frequency 
of the radiation (not the frequency of precession this time), and r is the length of time over 



which the neutrinos are released. (For burst radiation this would be a factor of 




4 Simple Dipole 

The 5-function is of course an idealizaton. The excess of neutrino energy in one direction is 
caused by the magnetic field and is expected to be of the form 

/(f]',t) = ^ (l + 3ecos^) . (39) 

This distribution is chosen so that it has only a uniform and a dipole contribution. The 
numerical factors are chosen so that / fdVL' = 1 and / f cos 6dfl' = e. Here e is again the 
fractional excess of neutrino momentum in one direction compared to the other, and 6 is the 
angle to the magnetic axis. 

If the magnetic axis points in the {6, (p) direction, and we are looking in the {6', cj)') 
direction, then it is easy enough to find 

cos 6 = n{9, 0) ■ n{9', 0') = sin 6 sin 6''(cos cos 0' + sin sin 0') + cos 9 cos 9' . (40) 

If we now insert equations p9|l and (P0|) into equations (j27j) and (j28|l each term of the 
integral over 0' will take one of the following three forms: 

/'27r /'27r /'27r 

/ cos20'(i0'= / COS 20' sin 0'(i0' = / cos20'cos0'(i0' = 0. (41) 
JO Jo Jo 

All three are zero, so the result for 9^^,qtt and 9^^,ott is identically 0. This is a reflection 
of the fact that there is no gravitational dipole radiation. 

The first term that will produce noticable radiation is the quadrupole term which we will 
write in the form 

f(n',t) = ^cos'9. (42) 

This is actually not a true quadrupole as it has a monopole contribution, but that is irrelevant 
as monopoles will not create GWs. The factor of An is to remain consistent with equation 
fl39|) . The factor of 7 represents the strength of this term. It is left general, but we would 
expect that if the dipole term is suppressed by the factor e, the quadrupole term should be 
suppressed by two factors of e so that 7 ~ e^. 
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If equation (ji^ is inserted into (|77|) and (j2HI) then we get 

r^oTT(x,t) = 2r-V(t-r)-sin2^cos20, 

6 

r^oTT(x,t) = 2r-V(t-r)Jsin2^sin20. 



(43) 



(44) 



This is nearly the same result as the 5-function case except for the numerical factor and a 
sin^^ instead of 1 + cos^. This small change results in suppressed radiation whenever the 
predominiant neutrino direction is alligned with the observer's line of sight. (Because it is 
a quadrupole term, there is no distinction between aligned and ant i- aligned.) Even more 
interesting, if we allow the axis of rotation to be different from the z-axis as we did in the 
case of the ^-function, we find that we only have constant terms and terms with angular 
frequency uj and 2uj. All of the higher frequency terms are absent. Unfortunately, it seems 
at first glance that 7 ~ so this radiation is likely to be suppressed by two more orders of 
magnitude and may be unobservable. 

5 Realistic quadrupole from an off-centered distribu- 
tion 

The best chance of observing such gravitational waves is from an object with a large 
quadrupole moment. The pulsar kick mechanism based on neutrino conversions |21 13 E] 
predicts an asymmetric neutrino emission whose strength and direction are determined by 
the magnetic field inside the neutron star. The magnetic fields deep inside the neutron stars 
are believed to grow during the first seconds after the supernova collapse {a — VL dynamo 
effect ^0]) because of convection. The cooling of the outer regions of the star creates the 
temperature and entropy gradient that causes convection in these outer regions. It is likely, 
therefore, that the strongest magnetic field develops away from the central region and that 
the overall magnetic field has a large non-dipole component. Since the neutrino emission is 
affected by this magnetic field, the asymmetric part of the neutrinos can form an off-centered 
beam. This, in turn creates a source for GW with a large quadrupole moment. 

Here we consider two cases, which correspond to resonant and off-resonant conversions 
of neutrinos in Refs. jS] and jjj, respectively. 

5.1 Resonant conversions 

Let us consider a -B-field of the form 




(45) 
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Here p is the distance from some magnetic axis, Bq is the magnetic field at the center of 
the star and po is the distance from the magnetic axis to the center of the star. This is the 
simplest model that has no singularity in the current at the magnetic axis and represents 
charged particles circling with a constant angular velocity so that the current increases 
linearly with p. This model is based on cylindrical symmetry so it is a simplification for our 
case, but it should be sufficient for an order of magnitude estimate. 

According to Ref. j21El5 the radius at which neutrinos escape can be described as 

r = ro + 5r, (46) 

The only change from the reference is that I'll allow B to vary with distance from the 
magnetic axis. For simplicity I'll place the magnetic axis in the x — z plane. Then we can 
find 5 as a function of 6 and 0: 

5r = 6q cos 9[1 — 2x sin 9 cos cp + x'^ sin^ 9] ; (48) 

6o is the same as 6 in the reference (except with the assumed constant B-field replaced by 
the B field at the origin.) The value x is simply the ratio of the unperturbed escape radius 
to the distance between the magnetic axis and the center of the star x = —. Both distances 
can be a significant fraction of the star's overall size and are not correlated so x is most likely 
of order 1 and can be greater or less than 1. 

The first term is the one discussed in Ref. |21 Ej it is responsible for the momentum kick 
and as demonstrated causes no gravitational radiation. The last term is mainly i = 3 and 
also doesn't contribute to the gravity waves. The middle term however, is a quadrupole term 
and does create gravity waves. 

Its worth noting that while this is actually a measure of how much the neutrinosphere 
deviates from the sphere, it is assumed that the neutrinos, responding to their environment, 
will have an energy spectrum that mimics the radius distribution. 

f(Q', t) = — (l + 3e(l + -x^)"^ cos mi - 2x sin ^ cos + x^ sin^ 9]] . (49) 
477 V 5 / 

This distribution has been chosen to have a large monopole term, and a perturbing term of 
the same form as equation (jlHjl . Like equation it has been normalized so that / fdfl' = 1 
and / / cos 9dQ' = e where e represents the fractional momentum excess in one direction. 

The next task is to find cos 6^ sin ^ cos in terms of (6*, 6*', 0, 0'). As I already discussed, 
cos^ is the projection of the r' {n') vector (the integration direction) onto the f (n)vector 
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(the direction of the magnetic field axis). Essentially, the radial direction is the z direction 
for the tilde variables, and cos 9 is the projection onto this axis. 

By contrast sin 6' cos is the projection onto the x-axis for the tilde variables. Since 
f is the z-direction, the x-direction has to be in the 6', plane, but it can be any linear 
combination of these. So I'll let it be 



X 



cos a9 + sin acj). 
In terms of the standard (x, ?/, z) directions we can write 



(50) 



e 



cos 6 cos ( 
cos 9 sin ( 
— sin 9 




(51) 
(52) 

(53) 



Which means 

sin 9 cos d) = cos a (cos 9 sin 9' cos 



— (f)') — sin cos 0' 

cos 6* = sin 6* sin 9' cos (0 — 0') — 
Substituting this into equation HHl and in turn into the GW equations (PTI andOH j) we get 



) + sin a (sin 6*' sin (0' — 0)) , (54) 
cos ^ cos (55) 



^-•,0TT (x,t) 

0-^OTT(x,t) 



2ea{t 



r (l + 
2ea{t 



2^2 
5-^ 



sin 9 



cos 9 cos a cos 20 — sin a sin 



cos 9 cos a sin 



+ sin a cos 20 / 

' cos (20 + 6) 



— , ^ sin 6'v cos^ 6* cos^ a + sin^ q; ( ^^"^^^T ,' \\ 1 • (56) 

r(l + |x2) ^ V sm(20 + 5) j ^ ^ 

Here tan b = tan a sec 9. 

Clearly there is radiation and there is no extra suppression compared to the dipole term. 
There is now a wider variety of radiation producing motions. Not only can we allow 9 and 
to change based on rotation about some axis, but a can also change. However, if we assume 
a simple rotation of angular frequency uj then a will be linked to the rotation angle. There 
can be a phase difference between them a = 0o + e, but e goes to if the magnetic and 
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rotation axes are coplanar or 7r/2 when maximally skewed. There is now an extra factor of 
the rotation angle which adds significant radiation at a frequency of 3uj as well as u and 2uj. 

The maximally skewed case is to be expected if the neutrino kick is also driving the 
neutron star's rotation as in ^1]. This could also suppress the translational kick and make 
the gravity waves larger than expected relative to the proper motion. 

5.2 Off-resonant neutrino conversions 

If the active-to-sterile neutrino conversions occur off-resonance, the sterile neutrinos are emit- 
ted from the entire core of the star Such emission is only possible after the electroweak 
matter potential is driven to zero by neutrino oscillations ^21- Therefore, in this case, the 
sterile emission can begin several seconds after the on set of the supernova. More impor- 
tantly, this emission may last for a period of time that could be as long as 10 s, but can also 
be much shorter. This depends on the neutrino mass and mixing parameters. 

The emission of neutrinos in the solid angle dil is ^ = Nq(1 + ecos O). The dependence 
on the magnetic field is more complicated than in previous models. However, we can get a 
reasonable order of magnitude estimate by parametrizing the emission in the same way as 
(j48|) . Specifically the gravity wave generating term will take the form 



This must be integrated over all radii where neutrinos are emitted. Assuming they are 
uniformally emitted throughout the core (i.e. no explicit dependence on radius) the final 
form will be 



Except for the extra factor of 3/4 the result is the same as for the resonant case and will 
yield a similar magnitude signal. 



We are now ready to consider the magnitude of the radiation. In a typical supernova, 
neutrinos carry away about 10^^ ergs over a period of several seconds. So a is approximately 
the average luminosity a = L = 10^^(10s/r) erg/s where r is the total time of neutrino 
emmision. The magnitude of the gravity waves is then 



— 2NoeoX sin 6 cos 6 cos (j). 



(57) 



NoeoXmax siu 6 COS 6. COS 



(58) 



6 Signal Magnitude 




(59) 
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F{9i,t) is the function of angles appearing in equation (f^ . 9i merely represents the collec- 
tion of angles. 

To be compared with experiment, this result needs to be integrated with respect to time 
and Fourier transformed. In the Fourier transform all three major components [u, 2u), and 
3a;) have amplitudes of order 1, and integrating is equivalent to dividing by the angular 
frequency. 

The units must be changed also. The gaussian noise must be Fourier transformed in a 
way that gives it units of B.z~^^^. To adjust for this we must divide by the square root of the 
frequency for a burst source, or multiply by the square root of r for a periodic source jTSj. 
Since neutron stars have typical frequencies of several kHz, and the burst is expected to last 
a few seconds, we should generally treat this as a periodic source. However, in the case of 
lower frequency signals detectable by LISA, the burst description may be more appropriate. 

The final result is 



for a burst signal. In both cases / represents the frequency of the radiation, not the star 
rotation. (Thus, / will be 1, 2, or 3 times the frequency of star rotation.) 

LIGO is most sensitive in the 10 - 1000 Hz range, neutron stars rotating at this rate 
are likely to undergo several rotations in the neutrino ejection time and therefore equation 
()60|) is most appropriate in describing LIGO's sensitivity to gravity waves figure ^ LISA by 
contrast is most sensitive in the .001 - .1 Hz range. Neutron stars rotating at this rate are 
likely to undergo only part of a rotation, or at most a few rotations in the neutrino ejection 
time. Thus, equation is more appropriate in describing its sensitivity figure |21 

From the figures it appears that LISA is more likely to see this GW signal, being sensitive 
to supernovas as far as several Mpc. Unfortunately, LISA's frequency range is very low and 
it would only see GWs from very slowly rotating neutrino beams. 




(60) 



for periodic signals, or 




(61) 



7 Gravitational Energy 




(62) 
(63) 
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10 20 50 100 200 500 1000 

frequency f, Hz 

Figure 1: Gravity wave intensity from a pulsar kick caused by neutrino oscilations. The 
LIGO sensitivity is shown for reference. For near supernovas with the appropriate rotation 
frequency, the signal should be observable. Please note this shows the magnitude at a given 
frequency. The actual signal is not continuous, but would be concentrated at three closely 
spaced frequencies, /, 2/, and 3/. 
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log f (Hz) 



Figure 2: Gravity wave intensity from a pulsar kick caused by neutrino oscilations. The 
LISA sensitivity is shown for reference. Note that for near supernovas with the appropriate 
rotation frequency, the signal can be observable. Please note this shows the magnitude at 
a given frequency. The actual signal is not continuous, but would be concentrated at three 
closely spaced frequencies, /, 2/, and 3/. 
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where 9^^ ,ott can be given by any of the several equations for this quantity as appropriate. 
Inserting equation (|59p this means that the energy carried away by gravity waves will be 
about .1 erg / cm^ or a total of about 10^^ erg. 

By comparison, Miiller and Janka [T3| find a GW energy oi Eq^^ ~ [10^^° — IO^^^JMqC^ ~ 
[10^^ — lO^^jerg, for a luminosity of ~ lO^^erg/s. Cuesta ^ gets a value of -Eqw ~ 
IO'^MqC^ ~ [10^°]erg with a luminosity of ~ lO^^erg/s. (The neutrino luminosity is 
important because gravitational luminosity scales with and the total energy scales like 
LIt ~ Ljy.) Thus, Miiller and Janka's estimate is of roughly the same size as mine, Cuesta's 
is about 3 orders of magnitude higher. 

Furthermore, the GW signal in this paper can be distinguished from other GWs caused 
by neutrino emission from a neutron star, by its strong periodic nature. This signal will be 
highly concentrated at three main frequencies /, 2/, and 3/ (see sectional). By contrast, 
other signals have their power spread out over a large, fairly continuous range of frequencies. 
In certain ideal situations individual cycles of the radiation may be recognizable. This should 
allow the periodic signal to be distinguished from other processes including those that contain 
more total energy. 



8 Conclusions 

We have shown that, in the event of a nearby supernova, gravitational waves caused by an 
asymmetric neutrino emission may be strong enough to be detected by LIGO and LISA. We 
have assumed that the asymmetry in outgoing neutrinos is as large as 1%, as suggested by 
the observed pulsar velocities. As can be seen from the figures, quickly rotating neutron stars 
may be seen by LIGO advanced, and more slowly rotating ones by LISA. This signal also 
has the distinguishing feature of having large radiation components at /, 2/, and 3/. This 
could be used for recognizing gravity waves from rotating neutrino beams. This mechanism 
could also be recognized by its short time duration and possible coincidence with neutrino 
and optical signals from the supernova. 

If such a signal is detected and identified, it would have profound implications for one's 
understanding of the supernova, as well as neutrino physics. In particular, LIGO and LISA 
may be able to test the existence of a keV sterile neutrino, which has been suggested as the 
explanation of both the pulsar kicks and the cosmological dark matter. 
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